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Abstract 

We investigate the nonleptonic decay of charmed meson into two pseudoscalar mesons 
using the vector-dominance model, and compare the results with those obtained from the 
factorization model. In particular, we discuss the role of the annihilation diagrams in the 
two models. 
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Much of what we understand about nonleptonic decays of strange particles comes 
from current algebra and the soft pion technique. Unfortunately, this technique cannot 
be applied to the decay of heavy mesons carrying the c or the b quark, since the emitted 
pions are generally not soft. In this case a completely different method has been developed 
based on the assumption that the hadronic matrix element of currents factorizes 1 . 

Historically, the nonleptonic decays of strange particles have also been discussed in 
terms of a dymanical model 2 based on the idea of vector dominance. In particular, it is 
well-known that for the K — > 2n decays this model provides a satisfactory description. 
It would be of interest to see how vector dominance fares in describing also the heavy 
meson decays. In this paper we deal with this question and, as an extension of K — > 27r, 
discuss the decay of D and Ds mesons into two pseudoscalars. Indeed such an analysis was 
attempted several years ago 3 . The present work obviously benefits from the availability of 
better data. More importantly, however, we undertake a detailed comparison of the vector 
dominance model and the factorization model. In particular, we discuss the role of the 
so-called annihilation amplitudes in the two models. 

For nonleptonic decay of charm, the effective weak Hamiltonian in the current-current 
form may be taken to be 4 



and ai, 02 are real coefficients which we treat as phenomenological parameters. The primed 
quark fields are related to the unprimed ones by the usual Cabibbo-Kobayashi-Maskawa 
(CKM) mixing matrix. We shall ignore the Penguin-type contributions. For the nonlep- 
tonic D and Ds decays into two mesons, the Hamiltonian (1) leads to two main classes of 
quark-model diagrams, the spectator and the annihilation diagrams shown in Figs. 1(a) 





where (q^q 01 )^ are color-singlet V-A currents 



(«,/?= 1,2. ..4) 
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and 1(b) respectively. It is well-known that the annihilation-diagram contribution in the 
quark model is helicity suppressed. 

The factorization model assumes that the matrix element for the decay D — > P\P 2 (Pi 
and Pi are light mesons) can be written in the factorized form 
< Pi( gi )P 2 (g 2 )|# w |D(g) > 

~^|[oi < P 2 (<z 2 )|M%|0 >< Palis' c)^\D(q) > (3) 

+ a 2 < P 1 (g 1 )|(5'd%l >< ^2(? 2 )|(«c) M |D(g) > ] 
The two terms on the right can easily be seen to correspond to the quark-model spectator 
diagrams of Fig. 1(a). An important feature of the factorization approximation is that the 
annihilation diagram can be neglected 1 . To see this, note that the annihilation diagram 
corresponds to the factorization 

^ a 2 < 0\(uc)^\D(q) X P 1 (q 1 )P 2 (q 2 )\(s'd') fl \0 > (4) 

if the decaying particle is a neutral D°. For the charge-carrying D + or Dg a similar 
factorization can be written down involving the a± term in the Hamiltonian. In these 
factorized forms, the charmed meson is connected to vacuum by an appropriate current. 
From Lorentz invariance, this matrix element is proportional to q^ = (qi + q2)n- Multiplied 
by the other matrix element in (4), we see that the annihilation contribution is then 
proportional to the matrix element of the divergence of the current (s'd')^ formed from 
light quarks, and is thus proportional to the masses of the light quarks. This is the 
analogue of helicity suppression in the quark-model annihilation diagram. By comparison, 
the spectator contribution (3) is proportional to quark masses involving the heavy charmed 
quark. 

There is another argument that shows an additional suppression of the annihila- 
tion contribution. To appreciate this, note that the structure of the matrix element 
< Pi(q , i)P 2 (q , 2 )|(s'(i / ) At |0 > in (4) can be written from Lorentz invariance in terms of 
form-factors to be evaluated at the momentum transfer (q± + q 2 ) 2 = q 2 = —m 2 D . In 
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the standard pole-dominated form, the form-factor is then expected to lead to a further 
suppression of the annihilation amplitude. By contrast, again, the form-factors appear- 
ing in the spectator amplitude (3) are to be evaluated at the low mementum transfer 
(q ~ <7i) 2 = 0.2 = ~ m 2 or (q ~ 92) 2 = q\ = and thus need not be suppressed. 

Consider now the vector-dominance model. Here we take the currents in Hw to be 
the hadronic currents given by the field-current identities (a, (3 = 1, 2 ... 4) 

(VJ% = y/2 gv(<l>J% 

(5) 

where (0 M )jg an d Pp are the field operators of the vector and the pseudoscalar mesons 
respectively, and gy and fp are the corresponding decay constants. The nonleptonic weak 
interaction (1) can then be represented by a two-meson vertex. For the parity violating 
decay D — > P1P2, the Cabibbo allowed and once-suppressed weak vertices are listed in 
Figs. 2(a) and 2(b) respectively. 

The Feynman diagrams for the decay D — > P1P2 are depicted in Fig. 3. It is easily 
seen that the diagrams in Fig. 3(a) and 3(b) are respectively the analogues of the spectator 
and the annihilation diagrams of the quark model. The annihilation amplitude in Fig. 3(b) 
is proportional to 

, S^ + q^q u /m v 

A a oc q^ 2- 2 — (91 - QV" (6) 

q + my 

In the pole term, if we set q 2 = —m 2 D , at first sight it seems to lead to a suppression in 
the annihilation amplitude. However it is trivial to see that the pole actually cancels, and 
one obtains 

. ml - ml 

A a oc 2 7 
m 2 v 

The annihilation amplitude does depend on the masses mi, m2 of the light mesons. Again, 
this is the analogue of the helicity suppression in the quark model. However, it is easy to 
see that V is a light vector meson in this case, so the annihilation amplitude is actually 
proportional to the ratio of light meson mass squares. By contrast, for the spectator 
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diagram of Fig. 3(a), the amplitude is proportional to 

A s oc (g + gi)^ a - ; 2 -?2^ = g 8) 

^2 + my 

This time since F is a charmed vector meson D* or D* s , we find that the spectator con- 
tribution while containing a piece proportional to the charmed meson mass term m 2 D , is 
actually determined by the ratio of heavy meson mass squares. Thus, there is no apriori 
reason why A a and A s may not be of the same order of magnitude, and we are not justified 
in neglecting the annihilation amplitude. 

We now use the vector dominance model to compute the decays D, Ds — > P\Pi-> 
taking into account contributions from both the Feynman diagrams in Fig. 3. We follow 
Bauer et al. 1 (BSW) and first determine the parameters ai, a 2 by confronting our model 
to the data on D — > Kn decays. As in BSW, we have to take the final state interaction 
into account, and we do this by considering only elastic scattering in the final Kit state. 
In terms of the isospin amplitudes, we have 

A(D° -> K~n+) = -L A 3/2 + ^1 A 1/2 

A(D° - K\°) = yjl A 3/2 - -L A 1/2 (9) 
A(D+ -> K°n+) = V3 A 3/2 

where 

A 7 = \A T \e iSl (10) 

is the amplitude in the isospin state I and 5i is the phase shift in that channel. Using the 
data from the recent particle properties data 5 booklet, we obtain 

\A 1/2 \ = 2.94 x 1(T 6 GeV 

\A 3/2 \ = 7.37 x 1(T 7 GeV (11) 
S1/2 - £3/2 = 93.4° 

We now use the vector dominance model to compute these amplitudes. Since we do 
not know all the strong interaction coupling constants, we use flavor SU(4) symmetry to 

5 



relate these to g p7T7T , which is known from the p — > 2n decay (g p777T ~ 4.0). The D 
amplitudes can then be written down from the Feynman diagrams of Fig. 3 to be 



Kn 



A(D° - 

A(D° 
A(D+ - 



R-7T+) =7 
+ K\°) 



m\-m 2 K m 

~ CL 2 gK*jD o 1" a l9D%l7i 
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(12) 



K°7i+) =7 



m l ~ m D , f m K - m D 
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where 
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i — V* s V ud 2g p7T7T 



V2 



(13) 



The V's are the CKM matrix elements and the couplings gy and fp in (12) are de- 
fined in (5). For the vector meson decay constant gv, we use the spectral function 
sum-rule 6 based on asymptotic SU(4), which predicts identical values of gyjmv for 
V = p, K*, D* 7 and D$. For the p-meson, we extract from p — > ££ decay, g p /m p ~ 0.15 
GeV. For the pseudoscalar meson decay constant fp, we take f n = 0.093 GeV, fx — 
1.2 fa = 0.112 GeV and choose fp, ~ fp> s ~ 0.136 GeV as determined 7 from the sum-rule 
technique in quantum chromodynamics. The isospin amplitudes in the vector dominance 
model can then be computed to be 

\A 1/2 \ = 1.07|(ai - 1.03a 2 )| x 10" 6 GeV 



\A 



3/2 1 



7.581 (ai + 1.34a 2 )| x 10" 7 GeV 



(14) 



Using the values (11) for these amplitudes obtained from the data, we find two solutions 
(only the relative sign of a\ and a 2 is important) 



I. a x = 1.98 



II. ai = 1.13 



a 2 = -0.75 



a 2 = -1.57 



(15) 
(16) 



It should be emphasized that the annihilation contribution has been included in this anal- 
ysis. We find however that in the D — > Ktt decays, this contribution is numerically small 
compared with the spectator contribution. 
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Isospin analyses similar to the one above can also be performed for the Cabibbo 
suppressed decay modes D KK and D — > tvtv. In the D — > nix decays, only D° — > tv + tv~ 
has a measured branching ratio at present. It is not hard to see that this experimental value 
is already sufficient to rule out the solution II. Also the surviving solution I is consistent 
with the data with negligible final state interaction. In the D — > KK decays, on the 
other hand, we find that the data require sizable final state interaction. In this case, while 
neither of the solutions can be ruled out, we find that I satisfies the data better. We do 
not present the details of these analyses here, but accept the conclusion that the present 
data favors the solution I over II. 

It is of interest to compare our solution to the one obtained by BSW in the factorization 
approach. Using the values (11) of the isospin amplitudes obtained from the recent D — > 
Kit data, we find that the BSW solution 8 is given by 

ai = 1.09 , a 2 = -0.48 (17) 

Our individual values of a± and a 2 in solution I are somewhat larger, although the ratio is 
not very different from the one in (17). 

We have calculated the various allowed and suppressed decays of the type D, Ds — > 
P1-P2, and have listed the results of the vector dominance model in Table 1. Also listed in 
the table are the results from the factorization model, and the experimental data, where 
available. 

In conclusion, we would like to emphasize that the vector dominance model affords a 
natural way of taking into account the annihilation contribution. While this contribution 
is small for the D, Ds — >■ PP decays, one would expect it to be larger in decays involving 
the heavier vector mesons as in D, Ds — > PV and D, Ds — > VV. In these cases, the 
vector dominance model has to be extended to include pseudoscalar and other meson poles. 
This work will be reported elsewhere. 

This work was supported in part by the U.S. Department of Energy Grant No. DE- 
FG-02-91ER40685. 
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Branching 


Vector-dominance 


Factorization 


Experiment 5 


Ratio 


Model 


Model 1 




BR(D° -> K~n+) 


4.9% 


5.0% 


(3.65 ± 0.21)% 


BR(D° -> Z°7r°) 


0.8% 


0.7% 


(2.1 ± 0.5)% 


BR(D+ -> 7fV+) 


2.7% 


2.7% 


(2.6 ± 0.4)% 


BR(L>+ -> iTif +) 


1.2% 


1.3% 


(2.8 ± 0.7)% 


BR(L>+ -> r77r+) 


3.7 % 


2.6% 


(1.5 ± 0.4)% 


BR(D° -> iT??) 


0.1 % 


0.3% 


<2.3% 










BR(D° -> K+K~) 


2.9 x 10" 3 


3.7 xlO" 3 


(4.1 ± 0.4) xlO" 3 


BR(D° -> K°F) 








(1.1 ± 0.4) xlO" 3 


BR(L>+ -> if +iT) 


7.4 x 10- 3 


9.6 xlO" 3 


(7.3 ± 1.8) xlO" 3 


BR(D° -> TT + TT") 


2.9 x 10" 3 


2.6 xlO" 3 


(1.63 ± 0.19) xlO" 3 


BR(D° -> 7r°7r°) 


0.8 x 10" 4 


2.5 xlO" 4 


<4.6 xlO" 3 


BR(L>+ -> 7r+7r°) 


1.4 x 10- 3 


1.1 xlO" 3 


<5.3 xlO" 3 


BR(L>+ -> K°7r+) 


4.6 x 10" 3 


2.3 xlO" 3 


< 6 xlO" 3 


BR(D| -> K+7r°) 


5.4 x 10" 4 


2.3 xlO" 4 





Table 1 
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Figure and Table Captions 

Fig. 1 Quark-model diagrams for the nonleptonic decay of D and Ds into two mesons. 

1(a) describe the spectator diagrams and 1(b) the annihilation diagrams. 
Fig. 2 Parity-violating weak vertices in vector-dominance model. 2(a) describe the Cabibbo- 

allowed vertices and 2(b) the once-suppressed ones. 
Fig. 3 Feynman diagrams for the decay D, Ds — > P\Pi in vector-dominance model. The 

weak vertex is represented by a black dot and the strong vertex by an open circle. 
Table 1 Branching ratios of the Cabibbo-allowed and once-suppressed decays of the type 

D, Ds — > P\Pi- The calculated values in the table ignore the final state interactions. 

For the vector-dominance model, we have used the values of a\ and a-i given in solution 

I. For the factorization model, we have used the BSW values given in (17). 
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